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Light scalar fields are expected to arise in theories of high energy physics (such as 
string theory), and find phenomenological motivations in dark energy, dark matter, 
or neutrino physics. However, the coupling of light scalar fields to ordinary (or dark) 
matter is strongly constrained from laboratory, solar system, and astrophysical tests 
of fifth force. One way to evade these constraints in dense environments is through 
the chameleon mechanism, where the field's mass steeply increases with ambient 
density. Consequently, the chameleonic force is only sourced by a thin shell near the 
surface of dense objects, which significantly reduces its magnitude. 

In this paper, we argue that thin-shell conditions are equivalent to "conducting" 
boundary conditions in electrostatics. As an application, we use the analogue of the 
method of images to calculate the back-reaction (or self- force) of an object around 
a spherical gravitational source. Using this method, we can explicitly compute the 
violation of equivalence principle in the outskirts of galactic haloes (assuming an 
NFW dark matter profile): Intermediate mass satellites can be slower than their 
larger /smaller counterparts by as much as 10% close to a thin shell. 
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I. INTRODUCTION 

Cosmological observations indicate that the universe is accelerating today with about 
75 percent of the energy density in dark energyi>^. If the dark energy is a scalar field 
then its mass needs to be very light indeed, but its couplings to ordinary matter must 
be suppressed to avoid fifth force constraints. Indeed fifth force experiments such as the 
Cassini satellite experiment put stringent bounds on the gravitational coupling of nearly 
massless scalar fields^i"-. The chameleon scenario posits that a scalar field with gravitational 
strength couplings to matter could generate the present day acceleration, but evade fifth 
force constraints^. The properties of the scalar field depend on the ambient density. In 
particular its mass is density dependent. Cosmologically the chameleon mass can be of 
order the Hubble constant, allowing the field to be rolling on cosmological time-scales^. 
However, in the solar system the chameleon can be sufficiently massive so that it evades fifth 
force constraints. Indeed, chameleon theories have a non-trivial way of evading empirical 
gravitational constraints via the existence of a thin-shell mechanism. For sufficiently large 
objects the chameleonic force is almost entirely due to a thin-shell of matter just below 
the surface of the object, with the matter in the core of the object giving a negligible 
contribution to the force. Although the original papers assumed the coupling to matter was 
of gravitational strength, it was later realized that, due to non-linear effects, the coupling 
could be much larger whilst still preserving the properties of the chameleon mechanisnt^. 

The chameleon mechanism has been tested in many different situations and constraints 
placed on the parameters of the theory^"^. Gravitational tests in the solar system^ and the 
full cosmological evolution have been studied^. The effect it has on structure formation on 
sub-galactic scales^ and on large scale structure formation have also been investigated^. 
More recently the chameleon coupling to photon fields has been investigated and con- 
straints have been placed on the mechanism from laboratory experiments^ and astrophysical 
measurements^. Laboratory experiments have been designed to look for the chameleonic 
force between objects^. These experiments are similar, though distinct, to Casimir force 
experiments and should either detect the chameleonic force or severely constrain it in the 
near future. 

On cosmological scales, the introduction of a chameleonic force can lead to violations 
of the equivalence principle, as objects with shallow potential wells (or small velocity dis- 
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persions) fall at a higher acceleration than the larger object that have a thin shell^. The 
chameleon effect on CDM large scale structure can now be explicitly seen in numerical 
simulations of f{R) gravity models^"— 

One effect that has not been studied is the chameleonic back-reaction on galactic scale 
objects. In chameleon theories large galactic haloes should have a thin shell. In particular, 
the Navarro-Frenk- White (or NFW) profile models the dark matter distribution in galactic 
haloes^i. In this paper, we investigate the back-reaction (or self-force) of satellites mov- 
ing in the halo outskirts of an NFW profile, which can lead to observable violations of the 
equivalence principle. To do this, we use an analogy between chameleon thin-shell condi- 
tions and electrostatics, which enables an application of the method of images to fifth-force 
calculations (This analogy was also recently noticed and used to argue the presence of 
'lightning rod' effects close to non-spherical thin shells). 

The plan of this paper is as follows: In the next section, we review the chameleon mecha- 
nism, explain how its field depends on the ambient density, and how the thin shell mechanism 
works in this context. We then compare the thin shell mechanism with electrostatics in Sec. 
mil showing that there is an analogy. This analogy enables us to use the method of images 
to compute the chameleonic back- reaction (or self- force). We compute the self- force cor- 
rections to the fifth force between a test object and a body with a thin shell before going 
on to consider the NFW profile and the chameleon force in Sec. IIVI In Sec. |Vl we show 
that the circular velocity of intermediate mass satellites in the outer halo can be reduced 
significantly if the back-reaction is taken into account. Finally, Sec. |Vl] contains a discussion 
of our results and concludes the paper. 



II. THE CHAMELEON MECHANISM 



Chameleon fields appear in scalar-tensor theories of gravity. The action of the chameleon 
field in the Einstein frame is 



S 



'-fR--{d<pf-V{<p) 



(1) 



where Mp = (SvrG) is the reduced Planck mass. Matter couples to both gravity and the 
scalar field according to 



fJ.U J 1 



(2) 
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where ip is a matter field and A the conformal factor relating the Jordon and Einstein frame 
metrics. Notice that the scalar field couples to all matter species, including baryons. This 
conformal coupling gives an extra contribution to the Klein-Gordon equation. 

= - (3) 

where = In the approximation that matter is well described by a pressureless, 

perfect fluid this becomes 

= V^ + a<i>PmA{(j)) (4) 

where pm is the matter density. 

An immediate consequence is that the dynamics of (p is governed by the effective potential, 
which depends explicitly on the density, 

Veff = V{<j))+p^A{<j)). (5) 

If the bare potential is of a runaway form, 

n<t>) = (6) 

where A is a parameter with dimensions of mass and the conformal factor increases with 0, 
for example as 

A(0) = exp(/30/Mp), (7) 

then the effective potential has a minimum that depends on the density. Thus the physical 
properties of the field depend on the ambient density. 

From the above we can see that the mass of the chameleon is an increasing function of the 
density such that it can be massive in dense environments, but nearly massless cosmologi- 
cally. Also, the non-linearity in the chameleon equation of motion gives rise to the thin-shell 
mechanism, which is discussed in the next section. The net effect of these properties is that 
chameleon fields evade detection via gravitational tests for fifth forces, but can play the role 
of dark energy. 

III. THIN-SHELLS, ELECTROSTATIC CONDUCTORS, AND THE 

CHAMELEON FIFTH FORCE 



It is well-known from electrostatics that to find the backreaction of a charge on electric 
potential close to a conducting sphere one may use the method of images. Here we consider 
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the analogous problem for the chameleon mechanism, in that we consider the backreaction 
of a test body close to a thin-shell body that is a "chameleon conductor" . The method of 
images will be effective insofar as the conductor approximation is valid, which in turn (as 
we argue) implies that the shell is sufficiently thin. 

Let us first summarize the chameleon properties close to a compact, spherical body with 
a thin shell. The chameleon field inside a body of radius Rc with a thin shell is^ 

< r < R,o\\ : ^ 0c, (8) 
^roU <r < Kc. (pint - iTjj- + — 7TT7 ^ V>c, (yj 



3Mp V 2 r J 2Mp 

where Mp = (SttG)"^/^, 0c is the value of the field inside the compact body of density pc 
and radius Rc, while Rroii is the radius where the field starts to move from its value inside 
the compact body. To make the above two approximations have assumed: 

Mp\V^\-^ (5pe^^l^"^' is negligible, (10) 

< 1 — > e^'*'!^'^-^ ^ 1. (11) 

We can write the field at the surface of the body, i.e. r = Rc, as 

<Ps - <P{Rc) = i^(^ + ^)- + 0c. (12) 



3Mp V 2 Rc J 2Mp 



If we are in the thin-shell regime, then 

ARc ^ Rc — RtoW 
Rc Rc 

and using a Taylor expansion, eq. ( !T2|) reduces to 



< 1 (13) 



Inserting the Newtonian potential 



into eq. (lT4j) we obtain 



1 Mc_PcRl f,^. 

- MdlR-c - ml ^^^^ 

p ^ p 



3(3Mp^M \ Rc J ' ^^^^ 
which provides a criterion for comparing the value of the field at the surface relative to 

its interior in terms of the thickness of the shell. Whilst the above has been derived using 
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approximations, it has been checked analytically and numerically in the original papers, both 
for /3 = 0(1)^ and in the strong coupling regime^. Using f fT3|) . we see that the field values 
are nearly the same, i.e. the field is continuous across the thin shell. Therefore, ignoring 
the outside chameleon mass, the problem becomes analogous to that of a conductor held at 
fixed potential in electrostatics, where (f) takes a fixed value inside, 0c, while satisfying the 
Laplace equation outside the object. As a result, the back-reaction problem can be treated 
analogously to that of a conducting sphere in electrostatics. 

In obtaining the total force between a spherical body with a thin-shell and a test object 
without a thin-shell, we have three types of forces to consider. 

• The net gravitational force from a density distribution p(r) of the thin-shell body 
acting on a test object of mass m. This is given by 

Fg = -^V</.G(r), (17) 

where 0g is the gravitational potential satisfying: 

^ f^. (18) 

• The chameleonic force acting on a test object of mass m coupled with strength (3 to 
a chameleon field 0(r) of a body with a thin-shell of mass M. Here the relevant force 
isl 

= -]^^V0(r). (19) 
Inserting the external solution for a thin-shell rigid body of radius Rc 

B Me"™"""" 

0exi(r) = + (t>oc. r> R, (20) 

and ignoring the exponential factor (since rriaoRc ^ 1) the chameleon force (or fifth 
force) can be written as 

mM , , 

F^ = -a^f (21) 

where a = (S'^/AirMp and r is the distance between the centers of the bodies. We have 
introduced M as the reduced mass of the thin-shell body with radius Rc and mass M 

M=^^M, (22) 
Rc 
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where 

where = Mc/SirMpRc is the Newtonian potential, <^oo minimizes the effective 
potential outside the thin-shell body and 0c minimizes the effective potential inside 
the thin-shell body which can be obtained as: 

where we have used (jS]) with n = 1 in the above and the density far away from the 
body with the thin-shell is poo, whilst the density inside the body is pc- 

• The back reaction force between the thin-shell body and test object. As noted above, 
this can be obtained using method of images, yielding 

1 mm' . . 

where m' is the image mass of the test body located at r' inside the thin-shell body 
of radius Re- 

, mRr , R^ /„„\ 

m' = r' = (26) 

Note that we have neglected the chameleon mass exterior to the thin-shell body, since 
this mass depends on the density of the environment and the density outside the body is 
assumed to be small (this will change in the next section, where we assume a diffuse NFW 
profile). Accordingly, the total force acting on the test body m from the thin-shell body M, 
including both the fifth force and back reaction, is 

Ft = Fg + F^ + 6F. (27) 

To summarize, the test object responds to not only the fifth force from the thin shell, 
but also to its own image refiected from the thin shell. While the primary fifth force is 
attractive, the back-reaction (or self-force) is repulsive, as the image has the opposite charge 
(i.e. negative mass). Moreover, the self- force will become comparable to the primary force 
as test object approaches the thin shell. 
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IV. NFW PROFILE AND CHAMELEON FORCE 

In this section we will obtain the chameleon force for a dark matter halo (similar to that 
of the Milky Way galaxy) acting on a test body by considering the NFW profile for the halo 
density^, which is of the form 



Pir) = ZTIT^ (28) 



where = lOkpc = 1.57 x 10^^ eV ^. The quantity ps is given by 

'Air 



300 kpc 

A-nr"^ p{r)dr = 200pc 







^(300 kpc)^ 



(29) 



where pc = SH^M^ is the citical density of the universe. Using the present Hubble constant 
Ho ~ 71 km/s/Mpc = 1.5 x 10~^^eV and the reduced Planck mass Mp = (STrG')"^/^ = 
2.43 X 10^7 eV we obtain 

Ps = 3.0 X 10"^ eV^ 

We also consider the potential as in Eq. ([6]). 

Next we divide the space into two distinct regions. 
I) The interior region. Here the following constraint 

^ « 1 (30) 

is satisfied, and so the chameleon equation becomes 

VV = + ^e^^/*'^^ (31) 

Neglecting and solving 

(32, 

where the exponential coefficient has been ignored since we have assumed /Scp/Mp <C 1, we 
find 

<Pintix) ^ 0, [x{l + xf] ^/^""^^^ , Xo< < X < Xoy (33) 

for the solution to eq. fl52]) . For convenience we have defined the dimensionless variable 

X = r/vg and 



l/(n+l) 

(34) 
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FIG. 1: Chameleon field (p versus distance x = r/vs, where the red solid line shows the analytic 
solution and the blue dashed one depicts the numeric solution. 

The quantities Xo< and Xo> respectively denote the lower and upper values of x that beyond 



which the constraint fl30|) is no longer valid. We can find these values by setting — 
and then solving numerically for the two real roots Xo< and Xo>. 
II) The exterior region. Here we impose the condition 

V.S < W 



0.01, 



which we take to be valid for x > Xc- The exterior solution can therefore be found by 
neglecting the potential term in the chameleon equation 

/3r2p(x) 



V (pextix) 



yielding 



with 



(pext{x) 



B\n{l+x) _C 

X X 



X > Xr 



B 



(35) 



(36) 



(37) 



where C and (poo are integration constants. 

The quantity 0oo can be interpreted as the minimum of chameleon field for x ^ Xc- 
However, we need some criterion for determining C and Xc- Recalling that at x = a;o> 
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the chameleon field starts to deviate from the interior constraint (130|) and that the exterior 
solution is valid for x > Xc, we find that if we demand 

^ « 1 (38) 

Xc 

then to a very good degree of approximation the interior solution is still valid inside the shell 
Xo> < X < Xc, and the shell can be considered as sufficiently thin. Hence to compute the 
integration constant C we match the interior and exterior solutions and their derivatives at 
X = Xc, which yields 

C = -5 In (1 + Xc) + (pooXc - <Ps [xc{l + Xcf] '^^"^'^ Xc, (39) 

where Xc is the smaller real root of the following equation 

^ ^ ' ^ 0, (n + A)xc + n + 2 ^ ' 

which can be solved numerically for given parameters n, P, A and (poo- 

For simplicity we henceforth set n = 1 and take /3 of order unity. Whilst we could 
work with different parameters and in the strong coupling regime, we would not expect our 
results to change qualitatively, as long as the thin-shell conditions are satisfied. Then the 
only parameters left to be determined are (poo and A. We determine these via two basic 
considerations. First, the thin-shell conditions f l5U]) and fl38p must be satisfied, where the 
former yields an upper bound for A of the order of 10 eV. Second, since the size of the halo 
is approximately 300 kpc then the values for (poo and A must imply that the smaller real root 
of eq. fHOjl be in a range between < -Rc < 300 kpc, where Rc = (10 kpc)xc. Numerically we 
find that the second consideration along with the upper bound on A imply the upper bound 
for (poo is of the order of 10~^Mp. Hence we choose (poo = 1-5 x 10~^Mp and A = 4.66 eV, 
which satisfies all considerations simultaneously, yielding Xc = 7.183^ . 

We now proceed to solve the whole chameleon equation numerically, using the information 
from our semi-analytic investigation thus far. Specifically, setting n = 1, (5 = 1 and A = 
4.66 eV yields 

M2x(l + a;)2 Mp"+2(/.(x)"+i 

while we have re-scaled (p/Mp — (p. We present our numerical and semi-analytical solutions 
in fig. (P . 
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FIG. 2: The effective mass of the chameleon field (red) and If (blue) versus distant x = r/rg. 
It is clear that for the internal region x < Xc = 7.183 we have nT-efj ^ If ^^S^l- 

It is evident that the numerical and semi- analytic solutions match with high accuracy, 
confirming the approximations we made for this choice of parameters. For example, one 
approximation neglected the mass of the chameleon field in the external region. Fig. ([2]) 
shows the effective mass of the chameleon, i.e. ^Ijf = d'^Veff{4>)/d<P'^, versus distant x. Here 
we see that the chameleon mass almost vanishes for x > x^, and that m^rf ^ |-^t^| for the 
interior region. The latter inequality is an equivalent expression for the interior constraint 
( 130|) : taking the derivate with respect to r from Eq. (!32|) gives fn^j-jr = —l3p'/(j)'Mp where (j) 
is the interior solution given by eq. ( l33l) . 

In the next section, we use our semi-analytic solutions ( l33l) and ( l36l) . to obtain phys- 
ical quantities such as the fifth force, its back-reaction, and their possible observational 
consequences. 

V. CIRCULAR VELOCITY FROM FIFTH FORCE WITH BACK REACTION 

In this section, as an application of the thin-shell framework developed above, we will 
obtain the circular velocity of orbiting satellite galaxies in the exterior region of the galactic 
haloes, caused by the gravitational force, the fifth force and the back-reaction. Since by the 
constraint ( l38l) we demanded the shell of the halo to be thin, we will use the method of 
images which taks into account the back-reaction effect of the satellites. 

Furthermore, for simplicity, we limit ourselves to those satellites in the galaxy that can 
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reasonably be considered as bodies without a thin-shell. For these satellites we have 

Qmp^N ' ^ ^ 

where (j)ext{x) is the external chameleon field (136|) calculated at the satellite location. Since 
the virial velocity of a satellite is v^^^ oc $Ar, the condition P2l) yields a constraint on the 
virial velocity 

4>ext{x) ^^^^ 



6(3Mp 

where we have neglected (pc, which is the minimum of the chameleon field inside the rigid 
body. For our choice of NFW halo and coupling parameters, this reduces to f^ir ^ 10-15 
km/s, close to the thin shell, which is a typical velocity dispersion for intermediate mass 
satellites of Milky Way. 

Inserting the NFW profile (l28i) for the density distribution of the halo in the eq. ([1] 
the net gravitational force (that is finite at the origin) is 



m 



B ^B\n{l + x) 



xil + X 



x^ 



(44) 



where B = psrl/2Mp. 

The fifth force acting on a test object of mass m caused by the thin-shell halo can be 
obtained from eqs. (fT9|) and ( l33l) for the interior region as 



/3m 



1 + 3x 



■'int 



or from 



/3m 



C 



x{l -|- a;)(l -|- n) 

51n(l + x) 



X < Xr 



B 



tC f 



(45) 



(46) 



x(l -|- x) x"^ 
for the exterior region. 

The force caused by the back-reaction effect is simply the radial force between the test 
object m and its image mass m', located at r' inside the thin-shell halo. This is given by 
eqs. ^ and ^ as 



5F 



m 



where we have used the dimensionless radial coordinate x = r /vg and Xc is the smaller real 
root in eq. (140|) . Finally, the total force is 



Fg + Fs + 6F, 



(48) 
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FIG. 3: Circular velocity versus distance from the center of the galactic halo. The green solid 
line is obtained by considering the net gravitational force while the red dashed one is obtained by 
considering the correction from including the fifth force. For the rest, the correction due to the 
back-reaction has been also taken into account. The mass of the test object increases from top to 
bottom: m = 10^ Mq (blue), m = W^Mq (cyan), m = 10^^ Mq (brown). 

which is the sum of the net gravitational force (1441) . the fifth force (I45p / fHB]) and the image 
force (iTj). 

From the radial force we may obtain the circular velocity v{r) = \J \F\r/m as a function 
of r, the distance from center of the halo. Including only the fifth force, it is easy to see that 
the circular velocity is the same for all the objects in the external region of the halo, i.e. it 
is independent of the mass of the test object. This changes once the back- reaction is taken 
into account, which is shown in Fig. ([3]). For example, we see that for an object with mass of 
the order lO^^M©, considering just gravitational force from NFW profile we will obtain the 
green line for the circular velocity. If we include the fifth force from the chameleon theory 
for the NFW profile we get the red dashed line, which has a negligible effect in the interior 
region, but manifestly changes the exterior behaviour. Including next the back-reaction we 
obtain the brown line, which clearly deviates from the dashed red one for a large satellite 
mass. For smaller masses the deviation is not as big: in fact, it is clear that the correction 
from the image force becomes less important as the mass of the satellite deceases. The 
region enclosed by the black triangle denotes the regime in which the back-reaction becomes 
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50% — 100% of the original fifth- force. Therefore, we cannot trust our approximations (which 
neglected the effect of satellite on the thin shell); the velocity should lie somewhere inside 
the triangle, which specifies the range of uncertainty in our prediction. This triangle region 
grows as the mass of the test body increases, a feature more vivid in the close-up. 

An interesting observational window into violations of equivalence principle on Galactic 
scales was introduced idM'^—, where it was argued that leading and trailing tidal streams 
of satellites of Milky Way would be asymmetric, if stars and dark matter (that dominates 
these galaxies) experienced different gravitational accelerations. Based on the symmetry 
of the tidal streams of the Sagittarius dwarf galaxy,— argue that the difference between 
these accelerations should be < 10%, which implies < 5% difference in circular velocities. 
Given that Wvir ~ 15 km/s, and M ~ 10^ Mq for Sagittarius dwarf^, we expect around 
3% difference in circular velocities (see Fig. [3]), which is just below the observational limit, 
assuming that Sagittarius orbit is just outside a chameleon thin shell. However, further 
improvements in these limits, and/or study of other tidal streams in the Milky Way halo 
could well provide a way to discover this effect in the outer Galaxy. 

VI. CONCLUSIONS 

In chameleon theories there is a fifth force, which is suppressed due to the chameleonic 
effect and the thin shell. In this work, for the first time, we have quantified the proper 
thin shell conditions in chameleon gravity for the dark matter haloes surrounding galaxies 
(such as Milky Way), assuming an NFW density profile. As a result we could obtain the 
chameleon force for such a profile, which adds to the net gravitational force. We showed 
that this additional fifth force manifestly changes the behaviour of the circular velocity for 
the test objects in the exterior region of the halo. 

Since chameleon gravity is non-linear, the addition of a test object can change the 
chameleon field, and thus lead to a back-reaction or self-force, an effect that has never 
been calculated explicitly. In this paper, we have calculated this back-reaction, or self-force, 
in chameleon theories on galactic scale objects. Our method used the analogy between 
gravitational objects with a thin shell and electrostatics, enabling us to use the methods of 
images to compute the back-reaction. We applied our method to the NFW profile of dark 
matter halos. 
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When we apply our results to the circular velocity of satellites in the surrounding region 
of the halo we have found that the back-reaction cannot be ignored. When we only included 
the fifth force we found that the circular velocity was the same for all objects in the exterior 
region surrounding the halo. This situation changes considerably upon including the back- 
reaction. Indeed, depending on the mass of the satellite, the back-reaction can be a large 
modification to the original result. 

Our results suggest that there could be a violation of the equivalence principle in the 
outskirts of galactic haloes. While current bounds, based on the observations of the tidal 
streams of the Sagittarius dwarf galaxy, are not yet sensitive to this effect, future surveys 
of kinematics in the Milky Way halo can dramatically strengthen these limits^. This opens 
up another channel for testing chameleon theories. More generally our results suggest that 
in some of the gravitational tests for chameleon theories the back-reaction should be taken 
into account when putting constraints on the parameters of the theory. This is currently in 
progress. 
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